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ACCURACY OF DISCRETE APPROXIMATION FOR INTEGRAL 
FUNCTIONALS OF MARKOV PROCESSES 

lU. V. GANYCHENKO, V. P. KNOPOVA, AND A. M. KULIK 


Abstract. The article is devoted to the estimation of the rate of convergence of integral functionals 
of a Markov process. Under the assumption that the given Markov process admits a transition 
probability density which is differentiable in t and the derivative has an integrable upper bound of 
a certain type, we derive the accuracy rates for strong and weak approximations of the functionals 
by Rlemannian sums. Some examples are provided. 

1. Introduction 

Let Xt, t > 0, be a Markov process with values in Consider an integral functional of the 
form 

( 1 ) iT{h)= rhixt)dt, 

Jo 

where h : R is a given measurable function. In this paper we investigate the accuracy of the 

approximation of /r(h) by the Riemannian sums 

71 1 

lT,n{h) = — X H^{kT)/n), n > 1. 

^ k=0 

The function h is assumed to be bounded, only; i.e., we do not impose any regularity assumptions 
on h. In particular, under this assumption the class of integral functionals which we investigate 
contains the class of occupation time type functionals (for which h = 1 a for a fixed A G !B(R'^)), 
which are of particular importance. 

Integral functionals arise naturally in a wide class of stochastic representation formulae and 
applied stochastic models. It is very typical that exact calculation of the respective probabilities 
and/or expectations is hardly possible, which naturally suggests the usage of the approximation 
methods. As an example of such a situation, we mention the so-called occupation time option 
[14], whose price is actually given by an expression similar to the Feynman-Kac formula. The 
exact calculation of the price is possible only in the particular case when the underlying process 
is a Levy process which is “spectrally negative” (i.e. does not have positive jumps, see [5]), and 
practically more realistic cases of general Levy processes, solutions to Levy driven SDE’s, etc. 
can be treated only numerically. To estimate the rate of convergence of the respective Monte-Carlo 
approximative methods, one needs to estimate the accuracy of various approximation steps involved 
in the algorithm. In this paper we focus on solving of such a problem for the discrete approximation 
of the integral functional of type (1). 

For diffusion processes, this problem was studied in [4] and recently in [11], by means of the 
methods involving the particular structural features of the process, e.g. the Malliavin calculus 
tools. On the other hand, in two recent papers [3], [13], an alternative method is developed, which 
exploits only the basic Markov structure of the process and the additive structure of the integral 
functional and its discrete approximations. One of the aims of this paper is to extend this method 
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for a wider class of Markov processes. To explain our goal more in details, let us formulate our 
principal assumption on the process X. 

X. The transition probability Pt{x,dy) of X admits a density pt{x,y) w.r.t. the Lebesgue 
measure on This density is differentiable w.r.t. t, and its derivative possesses the bound 


( 2 ) 


dtPt{x,y) 


< BT,xt ^qt,xiy), 


t < T, 


for some Bt,x > 1)/3 > and measurable function q, such that for each fixed t,x the 
function qt,x{') is a distribution density. 

In [3], [13], the condition similar to X was formulated with /3 = 1. Such a condition is verified for 
the particularly important classes of diffusion process and symmetric a-stable processes. However, 
in some natural cases one can expect to get (2) only with /? > 1. As the simplest and the most 
illustrative example one can take an a-stable process with drift: 


Xt — ct Zt, 

where c 7 ^ 0 and Z is a (e.g. symmetric) a-stable process. Then 

Pt{x,y) = t-^/^g(^^ 

where denotes the distribution density of Zi. Straightforward calculation shows that ( 2 ) holds 
true with /3 = max(l, 1/a), which is strictly greater than 1 when a < 1. Since the Levy noises are 
now used extensively in various applied models, the simple calculation made above shows that it 
is highly desirable to extend the results of [3] and [13], which deal with the “diffusive like” class of 
processes satisfying X with /3 = 1, to the more general case of X with arbitrary /3 > 1. 

Another aim of this paper is to develop the tools which would allow us to get the bounds of the 
form (2) for a wider class of solutions of Levy driven SDEs. One result of such a type is given in 
the recent preprint [10], with the process X being a solution of the SDE 

(3) dXt = b{Xt)dt + a{Xt-)dZt 

where Z is a symmetric a-stable process. The method used therein is a version of the parametrix 
method, and it is quite sensitive to the form of the Levy measure of the process Z on the entire 
Recently, apart from the stable noises, various types of “locally stable” noises are frequently 
used in applied models: tempered stable processes, damped stable processes, etc. Heuristically, 
for a “locally stable” process its “small jumps behavior” is the same as for the stable one, but 
the “large jumps behavior” of the former is drastically different from “tail behavior” of the Levy 
measure. Since (2) is genuinely related to “local behavior” of the process, one can expect that 
the results of [10] should have a natural extension to the case of “locally stable” Z. However, 
to make such a conjecture rigorous is a sophisticated problem; the main reason here is that the 
parametrix method treats the transition probability of a Levy process as the “zero approximation” 
for the unknown transition probability density pt{x,y), and hence any bound for pt{x,y), which 
one may expect to design within this method, is at least as complicated as respective bound for 
the process Z. On the other hand, there is an extensive literature on the estimates of transition 
probability densities for Levy processes (e.g. [1], [ 8 ], [9], [12], [15], [16], [17], [18], [19], [21]; this list 
is far from being complete), which shows that these densities inherit the structure of the densities 
of the corresponding Levy measures. In particular, in order to get the exact two-sided bounds 
for pt{x,y) one should impose quite non-trivial structural assumptions on the “tails” of the Levy 
measure even in a comparatively simple “locally stable” case. Motivated by this observation on one 
hand, and by the initial approximation problem which suggests the condition ( 2 ) on the other hand, 
we pose the following general question: Is it possible to give a “rough" upper bound, which would 
be the same for a large elass of transition probability densities of “loeally stable proeesses", without 
assuming complieated eonditions on the “tails" of their Levy measures? The answer is positive, and 
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it roughly says that one can get the bound (2), where at the left hand side we have the transition 
probability density of the SDE driven by a “locally stable” process, and at the right hand side we 
have a (properly shifted) transition probability density of an a-stable process. This bound is not 
necessarily precise: the power-type “tail” of the a-stable density might be essentially larger than the 
“tail” e.g. for exponentially tempered a-stable law. The gain is, however, that under a mild set of 
assumptions we obtain a uniform-in-class upper bound, which is clearly easy to use in applications. 
To keep the exposition reasonably compact, we treat this problem in a comparatively simple case of 
a one-dimensional SDE of the form (10), see below. The extension of these results to a more general 
multidimensional case is much more technical, and we postpone it to a separate publication. 

The structure of the paper is the following. In Section 2 we formulate and prove two our main 
results concerning the accuracy of the strong and weak approximations of an integral functional 
by Riemannian sums, provided that condition X is satisfied. In Section 3 we outline a version of 
the parametrix method, which makes it possible to obtain (2) for solutions to Levy driven SDEs 
without strong structural assumptions on the “tails” of the Levy measure of the noise. In Section 4 
an application for the price of an occupation time option is given. 


2. Accuracy of discrete approximation for integral functionals 

In this section we will prove two results. The first one concerns the “strong approximation rate”, 
i.e. the control on the Lp-distance between /t(^) and its approximation Ix^nih)- 

Theorem 2.1. Suppose that X holds. Then for any p > 0 


where 
(4) DT,g{n) = 


Ej 

= sup,^ \h{x)\, 

log n, 


p\ i/p 


lT{h)-lTAh) <CT,p||/i||(DT,/3(n))i/2, 


n 


I 1 

maxi 1,-^ 


n 


-1//3 


/3 = 1, 
/3>1, 


Ct,p = 


I (14p(p - l)BT,xy'^^T, p>2, 
\Ct, 2 = (28BT,xy/^T, pe(0,2). 


Remark 2.1. This theorem extends [3, Theorem 2.1], where it was assumed that /? = 1. 

The second result concerns the “weak approximation”, i.e. the control on the difference between 
the expectations of certain terms, which involve Ixih) together with its approximation lT,n{h). 

Theorem 2.2. Suppose that X holds. Then for any A: € N and any bounded function f we have 


(5) 


E,(/t(/i)) V(^t) - 


< 2^^‘^k‘^BT,xT^^^ 


'VWDtA^). 


Remark 2.2. This theorem extends [13, Theorem 1.1], where it was assumed that (5 = 1. In the 
proof below, we concentrate on the case /3 > 1. 

Using the Taylor expansion, one can obtain directly the following corollary of Theorem 2.2. 

Corollary 2.1. Suppose that X holds, and let ip be an analytic function defined in some neighbour¬ 
hood of 0. In addition, suppose that the constants D^p,R^p > 0 are such that 


<^(™)( 0 ) 


m\ 


< Dp 


-— , m > 0. 


Then for any bounded function f and a function h such that T\\h\\ < R^p, we have the following 
bound: 

^MlT{h))f{XT) - ^MlT,nmf{XT) < CT,X,h,Jf\\DTAA, 


CT,X,h,p = 2^^'^Dp,BT,X 




Ru 


1 + 


(-f) 


where 


-3 
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Before proceeding to the proof of Theorem 2.1, we give an auxiliary result this proof is based 
on. This result is, in fact, a weaker version of Theorem 2.2 with k = 1 and / = 1, but we give it 
separately to make the exposition more transparent. 

Proposition 2.1. Suppose that X holds. Then 

^xlrih) - Ea;/'r,n(/i) < 5BT,xT\\h\\DT,i3in). 

Proof. Let us introduce the notation used throughout the whole section: for t € [kT/n,{k + 
l)T/n),k > 0, we put r]n{t) = Cn{t) = fhat is, r]n{t) is the point of the partition 

{Tk/n, A: > 0} of the time axis, closest to t from the left, and Cn(A) is the point closest to t from 
the right, which is strictly larger than t. 

We have 

^xlrih) — ExlT,n{h) = [ 'E,x[h{Xs) — h{X^^(^g^)] ds 

Jo 


[ [ h{y)[Ps{x,y) - Pn„(s)ix,y)]dyds 

Jo Jr‘>' 


— Ml + M 2 , 


where 


r^n.jsX'/n p 

Mi= / h{y)\ps{x,y) - pr^^(^s)ix,y)]dyds, 

Jo JR'i 

M2= / h{y)\psix,y) - prj^(^s){x,y)]dyds, 

Jk^^^Tju 

for some 1 < < n, which will be chosen later. We estimate each term separately. 

For Ml we have 

rK.pT/n r j. „ 

|Mi| < ||/i|| / / [ps{x,y)+pr,^(^s}{x,y)]dyds = 2 \\h\\T-J^. 

Jo Jr<^ ^ 

Further, using (2), we get 

IM2I < ||/i|| / / \ps{x,y) - pr,^(^s){x,y)\dyds 

Jkn^^Tln Jm.^ 

< ||/i|| / / / \duPu{x,y)\dyduds 

Jkj^^^Tju Jr}Ti{s) 

pT PS P pT PS 

< BT,x\\h\\ / / / u~^qu,x{y)dyduds = BT,x\\h\\ / / u~^ duds 

JK rT/ti Jrin(s) Jr’’’ j k„_aTln JVn(s) 


(6) 


^^pT/n Jrjnis) 

|■{i+l)T/n ps 

= BT,x\\h\\ ^ 
i=k^ » 


/■{i+l)T/n rs 

JiT/n J iT 


iT I n 

p{i+l)T/n p{i+l)T/n 


u ^ duds 

u~^ dsdu 


= Bt^xWH ^ / 

i=K,p 

rp n—1 p[i-\-l)T/n n-i 

< —BT,x\\h\\ ^ / u~^ du =—BT,x\\h\\ j 

n JiT/n n A 


u 


du. 


A—u JiT/n 

^ — ^n,j3 


..pT/n 


Now we finalize the argument. 
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1) If /3 = 1, put = 1, n > 1. Then we get 

|Mi| <2||/i||rn-\ 

IM 2 I < —[ u~^ du = BT,xT\\h\\n~^ logn. 

n jT/n 

2) If /3 > 1, put + 1. Then 

|Mi| < 2\\h\\T- -- < 2\\h\\T - - -< A\\h\\Tn-^/^. 

To estimate M 2 observe that 


( 7 ) 


Therefore, 



rpl-P 


\ n ) “ /3- 1 


n 


1-1//3 


n 


1-/3 


< 




n 


1-1//3 


IM 2 I < —i?T,x||/i|| [ u ^ du < 

n Jk,pT/^ 


Bt,x 


T^-^\\h\\n- 


1 // 3 . 


□ 


Proof of Theorem 2.1. Since we can obtain the required bound for p < 2 from the bound with p = 2 
by the Holder inequality, we consider the case p>2 only. 

Define 

Jt,nih) := It{h) - It,n{h) = f Anis)ds, A„(s) := h{Xs) - 

Jo 

By definition, the function t 1 —)> Jt,n{h) is absolutely continuous. Then using the Newton-Leibnitz 
formula twice we get 


Therefore, 

where 


•7T,n(/l) ^ = p(p - 1) j Js,nW ^ ^n{s) ( / A„(t) dt ) ds. 


JT,n{h) "<p{p- 




Hi, 


P-2 

rCnis) 


l^n(s)| 

/ An{t)dt 

J S 

p-2 ( 


Js^n\h) 

An(s) 1 

/ An{t)dt\ 

J Cn{^) J 


ds, 


Let us estimate separately the expectations of H)p^p{h) and H^^p{h). By the Holder inequality. 
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Further, observe that for every s the variables 

An{s), \Js,n{h)r^An{s) 

are lF^^( 5 )-measurable; here and below > 0} denotes the natural hltration for the process X. 

Hence, 


/ p—2 

= Ea; ( J Js^n{h) A„(s)Ea 


'in (s) 


Anit) dt 


^Cn(s) 1 ^'5 


< Ea 


f Js,n{h)^ l^n('S)| Ex ( / An{t)dt 

Jo yJCnis) 


3-, 


(n(s) 


ds . 


By Proposition 2.1 and the Markov property of A, we have 


Ea 


'Cn(s) 


An(t) dt 


3^, 


Cn{s) 


Ex. 


Cn(s) 


L 


T-Cnis) 


An(t) dt 


< 5BT,xEDT,j3{n 


Therefore, using the Holder inequality, we get 


Ea;Ff|„p(/i) < 5Bt,xEDT, i3{n)\\h\\ ^E^, j Js,n{h) 


ds 


l-2/p 


< 10BTxT^^^^^DT,Jn)\\hf 


ds 


Ea 


1-2/p 




2/p 


Note that n ^ < Dx^Jn)., hence the above bounds for and Ea;Ffj,^p(/i) finally yield 

the estimate 


(8) Ea 


/ T 

JT,n(/i)f < 14p(p - l)BT,xT^^‘^'^DT,Jn)\\hf f E, ^ | 


ds 


l-2/p 


It can be easily seen that the above inequality also holds true if JT,n{E) in the left hand side is 
replaced by Jt,n{h). Taking the integral over t G [0,r], we get 


E. 


■nJXhfdt < 14p(„- (e. r 


ds 


1-2/p 


Because h is bounded, the left hand side expression in the above inequality is finite. Hence, 
resolving this inequality, we get 


E. 


: r JsAh) "ds < {lAp{p-l)BTxY/^T^^\DT,p{n)Y/^h\\J 
Jo 


which together with (8) gives the required statement. 


□ 


Proof of Theorem 2.2. Denote 

Sk,a,b ■= {(si, 'S 2 , •••, Sfc) G : a < Si < S 2 < ... < Sfc < 6}, /c G N, a, 6 G M. 
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We have 




(/T(h))" - fiXr) 


= k\ Et 


>SkA 




i=l 
k k 


= kl 


fn/i(2/*) j /(^) \^Psi-si.Ayi-i,yi)] PT-Skiyk,z)dzf[dyj n dsi 


>Sk, 0 ,T \i=l 


(9) 


'\i=l 

k 


j=l i=l 


-kl / YlHvi) f{z) YlPr,r,(si)-Msi.A){yi-i,yi) 

JSt,0,TJm'^+^ \t=l J ' 


\i=l 


k k 

Pt -r}n {sk)iVk,z)dzWdyj n dsi 

j=l i=l 

k 


= k\Y,j [ 

r=l dSk,o,T 


fe + 1 


n kiyi) I fi^)dsZ..,sk,Ti^^ yi,---,yk, z)dz JJ dyj dsi, 

/ ’ ’ ’ i=l i=l 


\i=l 


where the convention sq = 0 , Sk+i = T,yo = x, yk+i = ^ is used and the functions J^^\r = 1,... ,k 
are defined by the relations 

{yr-l,yr) - Prj„(sr)-Vn{sr-l)iyr-l,yrZ) ^11^^ i + l —Si {yi,yi+i)^ ■ 


^ ^Psr—S,— l[yr— 

Let US estimate the r-th term in the last line in (9). We have 




n kiVi) 1 yi,---, Vk, z)dz dyj dsi 


Sk, 0 ,T \~=1 


< 


'Sk,o,T 


j=l i=i 

k k 

kiiLsfe,T(^> yi,---,yk,z)\dzYldyj n dsi- 

j=i i=i 


Since the case /? = 1 was already treated in [13], for the rest of the proof we assume that /3 > 1. 
Consider two cases: a) Sr — Sr-i > kn^pT/n and b) Sr — Sr-i < kn,pT/n. 

In case a), using condition X and the Chapman-Kolmogorov equation, we derive 


/, 


j; 


(r) 


(Rti)fc + 1 

< Bt. 


si,...,Sk,T 


{x,yi, ■.. ,yk,z)\dzY{dyj 


1=1 


',X / Prin(sr-i)—r]n{so){X:yr-l) 


(“ — 1 


V ^qv,yr-i{yr)dv 


' f]n (^r ) (-^r—l) 

Since kn^i 3 > 2, then in case a) we have Sr — Sr-i ^ 2Tfn^ and hence 

T Sf — Sf \ 

Pni^Sfr) Tyjj(Sj._l) ^ Sr 'Sr—1 ^ ' 

n 2 


dyr-idyr- 
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Therefore, using the fact that qt,y{-) is the probability density for any t > 0 and y G we hnally 
get 


/ 

Jn 


J 


(r) 


(Kd)/=+l 


rSr—Sr-l 

<Bt,x v 

j — l J Sr-Sr-l-T/n 


~l^dv < 


Bt,xT ( Sr — Sr-l\-d 


In case b) we simply apply the Chapman-Kolmogorov equation: 

k 


L 


J 


(r) 


(Rd)fc +1 




{x,yi,... ,yk,z)\dzWdyj 


i=i 


< 

< 2. 


Prin{sr-i)-rin{so) (^) Ur-l) (j}sr-Sr-i iUr-li l/r) + Pri„{sr)-r]„{sr-i)iyr-li Ur)^ dyr—ldyr 


Therefore, summarizing the estimates obtained in cases a) and b) we get, using (7), the estimates 

k k 


■ ■ ,yk,z)\dz'[ldyj n dsi 

j=i i=i 


[ / 

JSk,o,T 

Bt,xT p ^.-1 

- n ^ l Jo 


(r-1) 


■{Sr Sr—l) 


(A: — r)! 


diS'f' — 


T rSr 


+ 2 


II 


Sr l (T - Sr) 


k—r 


< 


0 Jsr-K^pT/n{r -B)^- {k-r)\ 

(r - Srf-'' 


Bt,xT^, F 


-2f^ 


+ 


n 

2Tkn 


i 

Jo 


(r — 1)! {k — r)! 


dsr 


-dSr—\dSr 

rT 

kn,l3^/'^ 


^dl 


u ^du 


r'i 

- 

Jo 




n JO ~ !)• {k — r)\ 

- (/3-l)(A:-l)! +(A:-1)! 

^ 2^T^^^BT,xDT,y{n) 


{k-\)\ 

where in the fourth and the hfth lines we used that sj!zj < Sr~^- Taking into account that in (9) 
we have k terms, and the common multiplier k\, we finally arrive at (5). □ 


3. Condition X for solutions to Levy driven SDEs 
C onsider the SDE 

(10) dXt = b{Xt)dt + dZt, Xq = x, 

where Z is a real-valued Levy process. In [10] it was shown that if Zt is a symmetric a-stable 
process and b{-) is bounded and Lipschitz continuous, then the solution to equation (10) satisfies 
condition X with (3 = max(l, 1/a) (in fact, in [10] more general multidimensional SDEs of the form 
(3) are considered). In this section we outline the argument which makes it possible to extend the 
class of “Levy noises”. Namely, we will omit the requirement on Z to be symmetric, and relax the 
stability assumption, demanding Z to be “locally a-stable” in the sense we specify below. 

Recall that for a real-valued Levy process the characteristic function is of the form 

f> 0 , .^gIR, 
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where the characteristic exponent if: admits the Levy-Khinchin representation 

(11) V’(6 =-^«^ + + / (l - + iCwl{|„|<i})^((iu). 

^ Jn 

In what follows, we assume that = 0 and the Levy measure pi is of the form 

(12) pi{du) = + C-\u\~^~‘^ \ue{-ifl)du + m{u)du, 

with some C± > 0 and m{u) > 0 such that m{u) = 0 for |u| < 1, and 

(13) m{u) < |tt| > 1. 

On the interval [—1,1] the Levy measure p, given by (12) coincides with the Levy measure of a (non- 
symmetric) a-stable process. This is the reason for us to call Z a “locally a-stable” process: its 
“local behavior” near the origin is similar to those of the a-stable process. In that context condition 
(13) means that the “tails” of the Levy measure for p, are dominated by the “tails” the a-stable Levy 
measure. 

Let us impose three minor conventions, which will simplify the technicalities below. First, since 
we are mostly interested in the case /3 > 1, we assume that a < 1. Second, the latter assumption 
assures that the integral 

/ up{du) 

J{\u\<l} 

is well defined, and we assume that the constant a in ( 11 ) equals to this integral; that is, i/’ has the 
form 


m= [ (l-e*«“)Mdn). 
Ju 


Clearly, this does not restrict the generality because one can change the constant a by changing 
respectively the drift coefficient b{-) in (10). Finally, in order to avoid the usage of the Rademacher 
theorem (see [10, Lemma 7.4] for the case when b is just Lipschitz continuous), let us assume that 
b€C\R). 

In what follows we show how the parametrix construction developed in [10] can be modified 
to provide the representation and the bonds for the transition probability density pt{x,y) of the 
solution to (10) driven by the “locally stable” noise Z. 

Let us introduce some notation and give some preliminaries. We denote the space and the 
space-time convolutions respectively by 

{f*9){x,y)-= f{x,z)g{z,y)dz, 
jR<i 

if ® 9)tix,y) := / {ft-s* 9s)ix,y)ds = / / ft-s{x,z)gs{z,y) dzds. 

JQ do JRd 

Generically, the parametrix construction provides the representation of the required transition 
probability density in the form 

(14) Pt{x,y) =Pt{x,y) + f f p°t_s{x,z)^s{z,y)dzds, t > 0, x,yeR. 

Jo Jn 

Here p^{x, y) is a “zero approximation term” for the unknown pt{x, y), the function '^^(x, y) is given 
by the “convolution series” 


(15) 


'^tix,y) = '^^?'"{x,y), t>0, x,yGR, 


k=l 


the function 4>t(x,y) depends on the particular choice of p^{x,y), and equals 
(16) ^t{x,y) := {L^ - dt)pt{x,y), 


10 


lU. V. GANYCHENKO, V. P. KNOPOVA, AND A. M. KULIK 


where 

Lf{x) : = b{x)f'{x) + [ (/(x + u)- f{x))n{du), f G (E) 

Jr 

is the formal generator of the process X. The subscript x in above expressions means that the 
operator is applied with respect to the variable x. Note that to make the above construction 
feasible, one should properly choose the “zero approximation term” p^{x, y), so that the convolution 
series (15) converges and the space-time convolution in (14) is well dehned. To introduce in our 
setting such p^{x, y), and then to construct the bounds for the associated <ht(x, y) and its convolution 
powers, we need some more notation. 

Denote by the a-stable process with the Levy measure Pa,c±idu) = m*'"’^=‘=)(tt) du, 

^(«,C'±)(^) _ C+u~^~'^lu>Qdu + C'_(-ti)“^“"l„<o, 
and the characteristic exponent 

^(«,c±)(^)^ f (i_e*«“)^(“’^±)(du). 

Jr 

Note that since 

the process possesses the scaling property 

Law = Law , c > 0. 

Denote by the distribution density of By the scaling property we have 

Denote also by Z*-") the symmetric a-stable process; that is, the process of the form introduced 
above with (7+ = C- = 1. Let gl°‘^ be the respective distribution density and g^°J := g'f^\ 

Finally, denote by Xt{x) and 9t{y), respectively, the solutions to the ODEs 

(17) dxt = b{Xt)dt, xo = X, dOt = -b{et)dt, Oq = y. 

Note that these solutions exist, because b{-) is Lipschitz continuous. 

Now we are ready to formulate the main statement of this section. 

Theorem 3.1. Let 

(18) Ptix,y) ■= gt’^^\dt{y) - a;)- 

Then the convolution series (15) is well defined, and the formula (14) gives the representation of 
the transition probability density pt{x,y) of the process X. This density and its time derivative have 
the following upper bounds: 

(19) Pt{x,y) <C{gl‘)l\+gl''^){y-Xt{x)), t G (0,T], x,y G E, 

(20) dtpt{x,y) < Ct-^/'^{gll\ +gi°^){y - Xt{x)), t G (0,T], x,y G E, 

Consequently, the process X satisfies condition X with fi = 1/a and 

qt,x{y) = {gil\ + {y - xt{x)). 

Proof. First we evaluate $t(x,y). If it is not stated otherwise, we assume in all estimates obtained 
below that t G (0,T] for some T > 0, and x,y G E. Observe that G C^E). Indeed, this 

property easily follows from the Fourier inversion formula and the expression for the characteristic 
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function. It is known that — x) is the fundamental solution to dt — where 

denotes the generator of the process 


( 21 ) 

Since 

we have 


L(“-^±)/(x) = / (/(x + u) - f € C 2 ( 1 R). 

Jr 


dtPt{x,y) = 




fiA 


t2/« 

-(^(“>C'±)^(«.C'±)) 


w 

fl/o 


w=et{y)-x 


w 


( Ja,C+)y 
tVo^ > 


w 


w=0t{y)-x 


where in the last identity we used the scaling property of and the fact that is a 

homogeneous operator of order 1/a. Next, by the very definition of L and we get 

1 _^j^{a,C±)g{a)^ ( 


LxPt{x,y) = 


tl/o 


tl/o 


(a,c±)y 

^2/a ' 


W 


w=et(y)-x 


+ 


'Nl>i 


1 g{o^,c^)(y^ 






fl/a J ^ 1 /, 

Therefore, using the relation dtOt{y) = —b{6t{y)), we get 

( 22 ) 


w=et{y)-x 


X 


fl/c 


+ 


1 


(a,c±) ( Ot{y)-x-u \ _ ,^^c±) ( ^tiy) - X 


tllo 


tl/c 


{rn{u) — 


m 


{<^,c±) 


(it)) du 


4|>i 

=: ^l{x,y) + ^j{x,y). 

Further, we give the bounds for the absolute values of ^l{x,y), d>f(x,y), and ^t{x,y). In what 
follows, C denotes a generic constant, whose value might be different from place to place. One has 

(23) 5 ^"’‘"*^(x) < Cg^°‘\x), X G IR, 

(24) |(ff(“’^±^)'(a:)| < 0(1 + |x|)-iff(“)(x), X G R, 

(25) |(5(“’''±))"(x)| < 0(1 + |x|)-24^(x), X G E. 

Since the argument used in the proof of (23) - (25) is quite standard (see e.g. [10, Appendix A]), 
we omit the details. 

By (24) and the Lipschitz continuity of b{-) we have 

1^,■(...)! < (^) I < (^) ^ cp<»>(..(p) -.), 

To get the estimate for |<I>|(x,y)|, we first observe that 

\m{u) - m(“’^±)(u)|/|„|>i < Cg^'^\u), 


which implies 


^ C‘ 4 / 0 " ) A'(»)<ip + 


Otjy) - X 

tl/a 











































12 


lU. V. GANYCHENKO, V. P. KNOPOVA, AND A. M. KULIK 


Taking into account (23), we deduce that 

\^tix,y)\ < *^5"^ +g^^^){9t{y) - x) = C{g^^^{x) + g^^\x)) [6t{y) - x). 

Combining the estimates for <h^(x,y) and <h|(x,y), we get 

(26) \^t{x,y)\ < C{gll\+ gl'"^){et{y) - x). 

Our next step is to estimate the convolution powers of d>. It is shown in [10, Appendix B], that 
the kernel g[°‘\9t{y) — x) possess the following sub-convolution property: 

(27) f g^'^l{9t-siz) - x) 5 (^“) {9siy) - z)dz < Cg['^\9t{y) - x). 

JR 

By this property we get 


(28) 


[ gll\-s{^t-siz) - x)fi((“) {9s{y) - z)dz< Cg^^^{9t{y) - x), 

JR 

[ 9tl+i{9t-s{z) - x)gi‘^^{9siy) - z)dz < Cg[%{9tiy) - x) < Cg^^^{9t{y) - x), 
JR 


where in the last line we used that g^^ < Cgi, and therefore * 9^^ 

these estimates, we deduce in the same way as in [10, Section 3] that 


< Cg 


(a) 

t+2- 


Having 


(29) 


<hf(x,y)| < 


CoiCt) 


k-l 


k\ 


■{ 9 , 


(«) I „(“) 

1+1 + 9 t 


)[0,[v) 


x). 


Therefore, the series (15) converges absolutely for {t,x,y) G (0, 00 ) x IR, X H, and 


^t{x,y)\ < C{gll\ + gi‘^^){9t{y) - x). 


Applying once again the sub-convolution property (27), we see that the convolution ® T is well 
defined, and 

\{p° ®'it)^{x,y)\ < C{g[‘il\ + g[‘^^){9t{y) - x). 

Thus, the expression (14) for pt{x,y) is well defined for any {t,x,y) G (0, 00 ) x IR x R, and 

\pt{x,y)\ < C{gll\+g[‘^^){9tiy) - x). 

Finally, to get (19) we use the following inequalities, which were proved in [10, Appendix B]: 


(30) c\9tiy) - x| < \xtix) - y\ < C\9t{y) - x[. 

Since for any constant c > 0 we have g]^\x) >; gl'^\cx) for any t G (0,T], x, ?/ G R, this completes 
the proof of (19). 

Our final step is to use representation (14) in order to find the bounds for dtPt{x, y). Since (x, y) 
and $t(x,i/) are given explicitly, it is straightforward to show that these functions are differentiable 
with respect to t, and to check using (23) - (25) that 

(31) \dtpUx,y)\ < C'f^/"srJ"^(6»t(p) - x). 


(32) \dt^t{x,y)\ < Ct gi‘^^){9t{y) - x). 

To show that the convolution powers <h®^(x, y) are differentiable in t and to get the upper bounds, 
we use the following trick. The expression for ^f^^^^\x,y) can be re-organized as follows: 

<l>f<^’^+^\x,y)= [ [ <i>f^,{x,z)^siz,y)dzds 
Jo JR 

rtl2 


[ ^fXix,z)^,iz,y)dzds+ [ ^f\x,z)^t-siz,y)dzds 

Jo JR JO JR 


( 33 ) 
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If /c = 1, the first line in (33) does not allow us to differentiate ' {x,y), because the upper 
bound for dt^t-s{x, z) has a non-integrable singularity (t — at the vicinity of the point 

s = t (recall that a < 1). However, the identity given by the second line in (33) does not contain 
such singularities, and we can show using induction that for any k > 1 the function ^f^{x,y) is 
continuously differentiable in t, satisfies 

ft/2 n r-tl2 


®{k+l)i 


[•til p pt/z p 

\x,y)= / {dt^®’")t-six,z)<^siz,y)dzds+ / ^f’"{x, z){dt^)t-siz,y) dzds 

Jo Jr<^ Jo Jr'^ 

+ [ ^ft{x,z)^t/2{z,y)dz. 

jRd ' 

and possesses the bound 

(34) |a,4.f (i.i;)| < +g!“>)(9.(!<) -X). t > 1. 

Since the proof is completely analogous to the proof of [10, Lemma 7.3], we omit the details. 

From (34) we derive the following bound for the derivative of ^t{x,y): 

(35) \dt'i>f^{x,y)\ < + g['^^){9t{y) - x). 

Re-organizing representation (14) in the same way as (33), we get 

l‘t/2 p pt/2 p 

Pt{x,y) = Ptix,y) + / Pt_s{x,z)'i>siz,y)dzds + / p%x,z)^t-s{z,y) dz ds. 

Jo Jr'^- Jo JRd- 

Using the above representation of pt{x,y) together with (31) and (35), we derive the existence of 
the continuous derivative dtPt{x,y), which satisfies the inequality 

\dtptix,y)\ < Ct~^/°{g[l\+gi°'^){et{y) - x). 

Using estimates (30) in the same way as we did that in the proof of (19), we can change the argument 
^t(y) — a; in the right hand side of the above estimate to y — Xtix), which completes the proof of 
( 20 ). 


□ 


4. Application: the price of an occupation time option 


In this section, we consider an occupation time option (see [14]), with the price of the option 
depending on the time spent by an asset price process in a given set. Comparing to the standard 
barrier options, which are activated or cancelled when the asset price process hits some definite 
level (barrier), the payoff of the occupation time option depends on the time during which the asset 
price process stays below or above such a barrier. 

For instance, for the strike price A, the barrier level L and the knock-out rate p, the payoff of a 
down-and-out call occupation time option equals 

exp i^P hst<L}dt^ {St - A)+, 
and the price C{T) of the option is defined as 


C(r) = exp(-rr)A 


exp 


-P [ hst<L}dt] {St - K) 




where r is the risk-free interest rate (see [14]). 

Assume that the price of an asset S = {St, f > 0} is of the form 

St = 5oexp(Ai), 

where X is the Markov process studied in previous sections. Then the time spent by the process S 
in a set J C M equals to the time spent by X in the set J' = log J. 
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Let US approximate the price C(r) of our option by 


Cn{T) = exp{-rT)E 


n —1 


exp -pT/nJ2hs,r/n<L}dt {St - K)+ 


k=0 


Then using the results from the previous sections we can get the control on the accuracy of such an 
approximation. 

First we apply Theorem 2.1 and derive the strong approximation rate. 

Proposition 4.1. Suppose that the process X satisfies condition X, and assume that there exists 
A > 1 such that G'(A) := Eexp{XXT) = E{St)^ < +oo. 

Then 

Cn{T) - C{T)\ < exp{-rT)pG{X)^/^CT,x/iX-i){DTAn))^^\ 
with constants Ct^\/(x-i) and DtAA given by (4). 

Proof. The proof is a simple corollary of Theorem 2.1. Denote h{x) = p'i-x<\ogL- Then keeping the 
notation of Section 2 we get 

C{T) = Ee-^^^^\ST - K)+, CniT) = - K)+. 

By the Holder inequality with p = X and q = A/(A — 1), 


CniT) - C{T) 


< e 


-rT I 


E 


_ p—lT,n(h) 


X /( X - 1 )\ (•^-1)/'^ 


\e{St-K)X)"^'' 

Since for positive a and b we have |e““ — e~^\ < |o — 6|, then 

Cn{T) - C(r)| < e-^^G(A)i/^ {E\lT{h) - > 

and thus the required statement follows directly from Theorem 2.1 with p = A/(A — 1). 


□ 


We also can control the accuracy of the approximation using the weak rate bound from Theorem 
2.2. Observe that the bound given below is sharper than those obtained in the previous proposition 
precisely when A > 2. 

Proposition 4.2. Under the assumptions of Proposition f.l, we have 

Cn{T) - C{T) < 2f^^^+^max{BT,xpTHl + pT)exp{pT),G{X)}exp{-rT)DTAn), 


where 

DtAA ~ 

Proof. For some > 0 denote 

C^(r) = exp(-rr).F 

CAT) = exp{-rT)E 


(1 i/^)logn, /3 = 1? 


n 
max 


exp y-P {{St - K)+ A A^) 

exp (-pT/n'^I{s^j./„<L}dA {(St - K)+AN) 


k=0 


Then 


(36) 


Cn{T)-C{T) < cAt)-C^\T) + C(r)-C'"(T) + Cn{T)-CX{T) 


■^N 


N / 


^N/ 


We estimate each term separately. 
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Using Theorem 2.2 and the Taylor expansion of the exponent, we derive 

k\ 


OO ^ 

C^(r) - C^(r) < 2^^^BT,xNTeM-rT) ^ DT,^{n) 

k=l 


= 2f^^^BT,xpNT‘^{l + pT) exp(/3r - rT)DT,fi{n). 

For the last two terms we get 

C{T) - C'^(r) I + I C„(T) - C^(r) I < 2 exp(-rT).F [{St - K)+ - {St - K)+ A N] 

< 2exp{-rT)E[STl{ST>N}] 

STN^-\sr>N} 


= 2exp{—rT)E 

To complete the proof, put N = 


j^x-i 


. 2G(A) 

<j^exp{-rT). 


□ 
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